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We study both the background evolution and cosmological perturbations of anisotropic inflation- 
ary models supported by coupled scalar and vector fields. The models we study preserve the (7(1) 
gauge symmetry associated with the vector field, and therefore do not possess instabilities associated 
with longitudinal modes (which instead plague some recently proposed models of vector inflation 
and curvaton). We first intoduce a model in which the background anisotropy slowly decreases 
during infiation; we then confirm the stability of the background solution by studying the quadratic 
action for all the perturbations of the model. We then compute the spectrum of the hx gravitational 
wave polarization. The spectrum we find breaks statistical isotropy at the largest scales and reduces 
' to the standard nearly scale invariant form at small scales. We finally discuss the possible relevance 

I of our results to the large scale CMB anomalies. 
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I. INTRODUCTION 

The inflationary paradigm has become a widely accepted description of the early universe, which has been successful 



sptea 

in solving the classical problems of FRW cosmology 0,0 • Moreover, inflation provides a mechanism for generation 
I""!, of nearly scale invariant spectrum of perturbations that can seed the structure we observe today W\. Inflationary 
' era is characterized by a quasi de-Sitter expansion with a nearly constant Hubble rate. During inflation, quantum 
fluctuations are generated and amplified to scales above the Hubble scale (which specifies the horizon of casual 
processes) where they remain frozen until they re-enter the horizon after inflation has ended. A key assumption about 
the quantum fluctuations during the quasi de-Sitter stage is that they are in an adiabatic vacuum stage in the deep 
T"'. sub-horizon (early time) regime, which results in a nearly scale invariant spectrum. 
O ' Results of the WMAP experiment [1] are in overall agreement with the predictions of inflation. However, certain 
^ ; features of the full sky CMB maps seem to be anomalous in the standard picture. These anomalies include the low 
^ power in the quadrupole moment [sl], the alignment of the lowest multipole moments (known as the 'axis-of-evil') 
—i ' and an asymmetry in power between the northern and southern ecliptic hemispheres Q • The statistical significance 
of these anomalies has been extensively discussed in the literature and despite numerous efforts, an explanation due 
to a systematic effect or a foreground signal affecting the analysis is not forthcoming. These anomalies suggest a 
violation of statistical isotropy of cosmological fluctuations, which is at odds with standard mechanisms of generation 
' of perturbations in inflationary models. Therefore, the possibility of relating the anomalies with modifications in the 
, standard infiationary picture has been considered by numerous authors recently. 

' In the standard picture, aim coefficients of CMB temperature anisotropies satisfy < airnja* >~ Ci Sw 6mm' (so, 
\ the temperature fluctuations are said to be statistically isotropic) . It has been shown in [3| that an anisotropic 
I ' expansion that took place at the onset of standard single field inflation, leads to a nonvanishing correlation of the 
^\ aim coefficients for different /-modes. This violation of statistical isotropy might then be related to the alignment 
I of the lowest multipoles observed in the CMB data. Therefore, several authors considered inflationary models with 
^ ' an anisotropic stage. As the simplest possible modification, anisotropic initial conditions at the onset of single field 
• ^ , inflation have been considered in references [1, [^, [13, [IH ■ Isotropization of the universe takes place due to the fast 
' expansion supported by the inflaton field. Fluctuations at scales comparable to the horizon scale during the anisotropic 
^ , stage of inflation are sensitive to the background evolution and thus provide breaking of statistical isotropy at those 
- - • scales. If inflation had a limited duration, this could leave an imprint on the largest observable scales. On the other 
hand, fluctuations that leave the horizon after isotropization re-enter earlier than modes that left the horizon during 
the anisotropic stage, producing the standard nearly scale invariant spectrum at small scales ^ . Since in such models 
isotropy is reached very soon (within one Hubble time) the breaking of statistical isotropy can be visible if the total 
duration of inflation is tuned to only the minimal duration to solve the standard problems of FRW cosmology. 

The fine-tuning can be relaxed if the anisotropic stage is prolonged, due to the existence of other sources, for 
example vector fields. Recently, anisotropic expansion driven by vector fields have been studied by many authors. 
The anisotropic expansion is achieved when the vector field acquires a nonzero vacuum expectation value (VEV) along 
a spatial direction. There is a range of models which differ by the way the VEV is obtained. The first anisotropic 
inflationary model was considered long ago in reference [isj . where the U{\) symmetry of the vector field action is 



It was also shown in ,1^1 that the anisotropic stage could lead to a detectable gravitational wave signal. 
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broken by a potential term, which needs to have a negative curvature to support an anisotropic expansion. It has been 
shown in reference [3] that this model is plagued by a ghost instability, and its quantum theory is inconsistent ^. 
In another model, the ACW model [l^, the VEV is obtained by a lagrange multiplier field which fixes the norm 
of the vector field and breaks the U{1) gauge symmetry. The stability analysis of the ACW model by considering 
the most general perturbations of the background has been performed in references [l6l| and [13, [i3| • In the latter 
studies, it has been shown that the linearized per turbations diverge close to horizon crossing, indicating the instability 
of the background solution. Another model [l!l] considered a nonminimal coupling of the vector field to the scalar 
curvature which breaks the conformal invariance of the vector field action. The coupling with curvature allows for 
a slow-roll phase during the prolonged anisotropic expansion, and the universe isotropizes due to the existence of a 
massive scalar field, which is responsible for the overall isotropic expansion of the universe. Models of infiation with 
nonmininal coupling to curvature have also been considered for standard isotropic infiation [20| . Isotropy of space is 
realized if 3N mutually orthogonal vector fields have equal VEVs. Instead, for randomly oriented N vector fields, 
one expects a statistically isotropic background with order 1/Vn anisotropy. Perturbative calculations based around 
such a background configuration have been considered in references [U HI, 113 and in [s^] ■ The latter study is the 
only complete study linearized study of perturbations, taking into account all the physical degrees of freedom, and 
coupling between vector field and metric perturbations. It has been shown in [s^j that the equations of motion for 
the linearized perturbations become singular close to horizon crossing, leading to instabilities. In all of these models 
(including [l^, [l^ [131), instabilities are related to the existence of the longitudinal polarization of the vector field 
(which would otherwise be absent when the C/(l) symmetry is restored), which becomes a ghost close to horizon 
crossing (In the vector inflation model [13|, ghosts also appear in the deep UV regime). 

More recently, reference (2^ introduced an anisotropic inflation model driven simultaneously by a vector field and 
a massive scalar field. The vector field is massless, but it is coupled to the scalar field through its kinetic term. Such 
type of coupling preserves the C/(l) gauge invariance of the vector field, so that the dangerous longitudinal vector 
polarization is absent; moreover the conformal invariance of the vector field is broken. When the vector field has a 
nonvanishing VEV along a spatial direction, the model possess an anisotropically expanding attractor solution. The 
anisotropy in the attractor is initially small, but it increases towards the end of infiation; therefore this is a counter 
example to the cosmic no hair conjecture (See ^25,] for a different example) .The breaking of conformal invariance due 
to the scalar-vector coupling can also be used to generate magnetic fields from inflation, as discussed in reference [1^ 
and more recently in [13, [11| and [IJ. These models are expected to be free from ghost instabilities as long as 
the coupling to the scalar field remains positive. However to our knowledge, there has been no complete study of 
the stability of these models, which take into account all the physical degrees of freedom of the system (including 
gravity). A complete study of stability is therefore necessary, given the problems identified with other vector field 
models [3[I3,[3g. 

In this work, we will present a complete study of cosmological perturbations of models with scalar-vector coupling 
that lead to an anisotropically expanding Bianchi-I background solution. We will also consider generalizations of the 
original model introduced in [2J|, to include additional scalar fields, so that isotropization can be achieved before the 
end of inflation. This is required in order to obtain a scale invariant spectrum of perturbations at small scales, but the 
large scale spectrum is modified, which in turn can be related to the CMB anomalies. Our study has two main steps: 
one is to show that the model is free of ghost instabilities, and the second is to study the resulting phenomenology. 
We perform the first step explicitly in this paper. We develop the necessary tools to study the phenomenology of 
the model and as a simple exercise, we study only the spectrum oi hx gravitational wave polarization in this paper. 
The study of the full phenomenology (including the spectrum of the curvature perturbation) will be communicated 
elsewhere [Slf . 

We will follow the formalism developed in [1, [l3| to decompose and classify the perturbations. The generic back- 
ground metric we study is the Bianchi-I metric with a residual two dimensional isotropy, given by 

ds^ = -dt^ + aitf dx^ + h{tf (dy^ + dz^) 

We exploit the two dimensional isotropy of the y — z plane to decompose the perturbations into two decoupled classes: 
2d scalar modes and 2c? vector modes (As discussed in [^ , there is no two dimensional transverse-traceless mode) . The 
linearized Einstein equations and the quadratic action for the two types of modes are decoupled and therefore we study 
them separately. We can deduce the number of degrees of freedom coming from each of the fields (gravity-|-vector 
field-fscalar field) by a simple counting, which does not depend on the decomposition chosen to classify them. The 
metric has 10 perturbations {5g^i^) to start with, 4 of which can be removed by coordinate transformations. Of the 



^ One may hope that this model has a well behaved UV completion. However, inflationary predictions are sensitive to this high energy 
regime, where the completion is needed. 
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remaining 6 modes, 4 are nondynamical, which can be best understood from the ADM formahsm [s^]. In the ADM 
formalism, the gravitational action is decomposed into the dynamical part containing the spatial metric hij and a part 
containing the lapse (N) and shift functions (iV,). The lapse and shift functions have no kinetic terms in the action, 
and they can be integrated out by solving their equations of motion This leaves only hij as dynamical modes, which 
have only two degrees of freedom. In summary, the metric perturbations have only 2 dynamical degrees of freedom 
(which are the gravitational wave polarizations in the standard case) and 4 nondyncamical modes. For the case of 
the vector field, out of 4 perturbations to start with ((5A^), one perturbation can be removed by the U{1) gauge 
transformation. Out of the remaining 3 perturbations 2 of them are dynamical and one mode is nondyamical (SAq). 
The scalar field perturbation introduces a single dynamical degree of freedom. Therefore, in total, the perturbations 
comprise of 5 dynamical modes, 2 of which are 2d vector and 3 are 2d scalar modes. When additional scalar fields 
are considered, each field introduces an extra scalar dynamical degree of freedom. 

We will insert the perturbative expansion of the metric, the vector field and the scalar field(s) into the starting 
action and expand it at the quadratic order. We will show that the actions for the 2d vector and 2d scalar modes are 
decoupled, so we study them separately. We will determine the linear combinations of p erturbations that canonically 
normalize the action, generalizing the computation of the standard isotropic case [33|. The study of the quadratic 
action is crucial for both showing that the model is consistent (free of ghosts) and also for determining the initial 
conditions for the perturbations. For modes that are inside the horizon, canonical combinations of perturbations can 
be quantized, and initial conditions can be set by the canonical commutation relations and by the requirement that 
the adiabatic vacuum state has minimal energy. We then proceed by numerically integrating the evolution equations 
for the canonical fields, starting from adiabatic initial conditions deep inside the horizon, until the end of inflation, 
which gives us the primordial power spectra. We do so for the 2d vector modes in this paper (which will be related 
to the power spectrum of the hx polarization of gravitational waves) and perform the study of the spectrum of 2d 
scalar modes in a separate publication. 

The paper is organized as follows: In section [TTl we summarize the anisotropic inflationary background solution 
obtained in reference [IJ] and generalize the model to possess extra scalar fields. This generalization leads to isotropiza- 
tion before the end of inflation. In section IIIIl we discuss the perturbations around the background configuration. 
Specifically, in subsection llll Al we review and discuss the decomposition of perturbations around the Bianchi-I back- 
ground solution, following the formulation of reference [9] . In subsection IIIIBl we discuss the generic properties of 
coupled vector-scalar models and develop tools that we will use in the following sections in order to find adiabatic 
initial conditions for such systems. 

In subsection IIII CI we discuss 2d vector perturbations around the background configuration. We compute the 
action for 2d vector modes and find the combinations of perturbations that canonically normalize the action. In 
section IIVI we compute the power spectra of the 2d vector modes numerically. This study results in the spectrum 
oi hx gravitational wave mode. We show that the spectrum has angular dependence (so breaks statistical isotropy) 
at large scales and reduces to the standard nearly scale invariant form at small scales. We also provide a fit to 
the numerically obtained spectrum. Finally in section |V] we provide a general discussion about the results we have 
obtained and their possible relation to observations. 

We also provide two extensive appendices: In appendix IVI Al we derive the relation between the 2c? decomposed 
modes and the standard longitudinal mode, which has been used to determine the power spectra. In appendix lVI Bl we 
perform the study of the 2d scalar modes. More precisely, we compute the quadratic action and determine the modes 
that canonically normalizes the action. The spectrum of scalar perturbations and the consequent phenomenological 
predictions (specifically, the < a;,„ „j/ > correlation) will be presented elsewhere. 



II. ANISOTROPIC INFLATIONARY EXPANSION DUE TO COUPLED VECTOR AND SCALAR 

FIELDS 



The models we discuss in this work are generically described by the following action 
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The original discussion of these coupled vector-scalar field theories in an anisotropic inflationary background was 
given in [2^, where there is only the single scalar field </> in the action. We assume that only one of the scalar fields 



The equations of motion derived from extremozing the gravitational action with respect to lapse and shift functions result in the 
momentum and hamiltonian constraints. 
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(or only one linear combination) enters in the kinetic function / for the vector field. We denote the remaining scalars 
with (f>a (a = 1, . . . , N). The metric ansatz is taken to be the homogeneous but anisotropic Bianchi-I metric given by 



=2a(t) 



dx^ + e2-(*) {dy^ + dz^) 



(2) 



In the above metric a measures the number of e-folds of average isotropic expansion of the universe and a measures 
the deviation from anisotropy. This metric also possesses a left over 2d isotropy in the y — z plane (which we will 
make use of when decomposing the perturbations). We can also identify ^ with ds^ — —dt^ + a? dx^ + \P {dy'^ + dz^), 
then we would have 



„q-2ct 
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(3) 



where Ha = a/ a, Hb = b/b, H is the average isotropic expansion rate, and h is the anisotropic expansion rate. An 
overdot denotes derivative with respect to the cosmic time t. An isotropic flat FRW universe is the limiting case 
when h 0. We will make use of both representations of the Bianchi-I metric in the following discussions. The field 
equations obtained from the action ([1]) are 
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where 



(5) 



and /'(0) = df /d4>. We assume that all the scalar fields are homogeneous so that (j) — (j){t) and = 4)a{t) for all a. The 
vector field is assumed to have a homogeneous vacuum expectation value along the x-direction so = (0, ^i(i), 0, 0) 
(a homogeneous Aq does not enter into the field equations, and can be set to zero by using the U (1) gauge invariance). 
With the metric and vector field ansatz, the last of ^ can be integrated to give 



^1 = / ^ (0) PA e 



-a — 4(7" 



(6) 



where pa is an integration constant. We use ([6|) and the metric ansatz ([2]) in the field equations ([4]) to obtain the 
evolution equations 



where we have defined 
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The values of the functions a and a are nonphysical; what is physical is their time derivatives. Therefore physics is 
the same under a constant shift a{t) a{t) + Uq and a{t) —+ a{t) + ctq (This corresponds to a rescaling of the spatial 
coordinates by two constant factors.). Clearly, the value of the constant pa is nonphysical, and it changes under a 
coordinate transformation. We can deduce its transformation properties by imposing invariance of Fa(3 F"^ which 
reduces to the condition 

Fall F""!^ ^ F^p F°'l^ under {a ^ a + ao, a ^ a + ao} ^ , pa ^ PAe^'^"'>+''°'> (9) 

Thus, PA is a physical parameter invariant under the constant coordinate rescaling. For this reason, it is the combi- 
nation which enters in the field equations ([7]). 

We will discuss inflationary solutions of ([7]), with appropriate choices of f{(p) for both single field (only </>) and two 
field models (a=l) respectively. 



A. Single Scalar Field Inflationary Background 

In this subsection, we review and discuss the background evolution of the original model [131 with a single scalar 
field coupled to the vector field. As described in [2J], we look for solutions that have a slow-roll regime and small 
anisotropy. Thus, we neglect &, ip in the first and a in the last of Q to get 
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+ T^f (<!>) ' 3d</)«-l^'(0)-pi/-3(0)/'(</)) 



2M2 



(10) 



Up to now, f{(f>) has been arbitrary. We now specify it so to obtain the desired anisotropic background solution. To 
do so, we first assume that the effect of the vector field is completely negligible, and the resulting equations can be 
solved approximately in the standard slow-roll approximation, given by 



a : 



1 



(11) 



The case cr = corresponds to the isotropic FRW background. A prolonged anisotropic inflation requires Ai to be 
nearly constant. From equation ^ we see that this can be achieved if we choose 
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For instance, when V (x (/)"■, this becomes / ~ e 



Furthermore, one can set / ~ e 



(where c is a constant) 



and for c > 1, the anisotropy is expected to grow. From now on, we assume that the functional form of /((/>) is given 
by 



f{4>) = Exp 



" 2c 


1 ""^fl 


Ml J 





(13) 



We now return back to equations (|10p and look for solutions that have small anisotropy, but the vector field contribu- 
tion is not totally negligible. We specify the scalar field potential as V{(j)) = m? (jP' /2 and set f = e'^'^ /^^^p fj-om now 
on. The vector field modifies the evolution of the scalar field, if the terms V'{(j)) and p\ f~^{4') f'{4>) are comparable. 
Namely, for our choice of the potential and /((/>) 
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■4''' /Ml 
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For such type of solutions, wc can compare the energy densities of the scalar and vector fields. The ratio, (when the 
scalar field is still in the slow-roll regime and the anisotropy is small) is given by 
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M2 
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(15) 



where, the second approximate equality is obtained from (|14[) . This ratio is much smaller than 1 during inflation, 
since 4)/Mp ^ 1. Thus, even when the dynamics of the scalar field is modified by the action of the vector field, as 



6 



long as the anisotropy is kept small enough, the energy density of the vector field can be neglected. Therefore, we 
can rewrite equations pH)) by neglecting the effect of the vector field only in the first equation as 
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These equations are integrated to give 
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where D is an integration constant. When the second term in the right hand side of the above equality is dominant 
we get 



M2 
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(18) 



which is the standard slow-roll relation. The second term in the right hand side of ([T7]) will eventually be subdominant 
for c > 1, since the universe expands and a grows Therefore, the first term eventually dominates and we have 



2 ~2 

c Pa 
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2 mMp 

3 c 



(19) 



which has an extra 1/c factor compared to the standard slow- roll solution (|18p . Using the third of ([7]), for small and 
slowly varying anisotropy we have 



3 d (7 
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Moreover, combining this result with p9|) and the first of (flB)) we find that the anisotropy obeys 
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a 



9 M2 d2 



< 1 



(20) 



(21) 



which is always satisfied (either in the region described by ([T8|) or (fT9|) '). Also, note that the anisotropy is slowly 
increasing with time, since (j) is decreasing towards the end of infiation. The slow-roll solutions described by is the 
standard isotropic infiationary attractor and the solutions described by (jl9p is the anisotropic inflationary attractor. 
We are interested in configurations that start close to the attractor solution p^ . For illustration, we solve the system 
d?]) numerically by setting the initial conditions on the attractor solution. Namely, we set initial conditions as 
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• m 
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Inserting 



M2/2(0) 

into the constraint equation (the first of ([7]), we solve for din to get 

c-1 (bh 
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/25(c-l)2 2(1 -I- 3c) 



M2 



9c4 



3c2 M2 Mi 



(22) 



(23) 



For c < 1, the initial anisotropy will decay exponentially fast so we do not consider this possibility here. The solutions with c < 1 are 
disconnected from the anisotropic attractor solutions c > 1, that is there is no smooth limit of c — > 1 which takes us back to the isotropic 
FRW solution. We will discuss the implications of this effect on the perturbations in the next section. 
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FIG. 1: The left panel shows the evolution of anisotropy during inflation and the right panel after inflation (during when the 
scalar field is oscillating). H = a and h = &. The black dashed curve represents the approximate slow-roll solution during 
inflation and the straight red curves the numerical solution. 



ttCe-folds) 




FIG. 2: The left panel shows the evolution of the average Hubble rate {H = a) and the right panel is the phase plot of the 
scalar fleld. Again, the black dashed lines are the slow-roll solutions. Time is measured in units of m. 



Equations and specifies the slow-roll initial conditions, which are set on the attractor solution. The only 
free parameters are (pin, am and am- The values of am and ai„ are nonphysical, so we can set ai„ — (Jin — for 
convenience (which is equivalent to setting ain — bin — !)• Thus, only (pin needs to be specified, which determines 
the amount of infiation. The number of e-folds of inflation is approximately given by a{tend) — N Ki c(/)f„/4Mp, 
so for c = 2 we set (pin = llMp which gives around 60 e-folds of infiation. We show in Figs. [H [2] the numerical 
evolution of the system. Note that the slow-roll and the numerical solutions are in good agreement with each other. 
The anisotropy in the single field model is increasing towards the end of infiation and decreases after inflation (as 
can be seen in the right panel of Fig. [T|) during when the scalar field has an effective equation of state of matter. We 
are however interested in cases where the anisotropic expansion takes place at the onset of inflation (for example as 
in [1, [13) El and [l^). The main reason is because we only want the largest observable scales to be modified from 
the perturbations that left the horizon during the anisotropic initial stage. This in turn might have some relevance 
to the alignment of the lowest multipoles. In the current case, all observable scales are modified (indeed smaller 
scales are modified even more since the anisotropy increases towards the end of inflation) and this is not compatible 
with observations. A simple possibility to overcome this difficulty is to introduce extra scalar fields as in the starting 
action ([1]) . We discuss the background evolution for the simplest possibility of a two-field modification of the original 
proposal of in the next subsection. 



B. Two Scalar Field Inflationary Background 

In this section we discuss the simplest multi- field case of the action ([T]) with a = 1. In this model, the extra scalar 
field (pi is not coupled to the vector field, and it is only relevant for the overall isotropic expansion of the universe. For 
the rest of the paper we will assume that V{(p) = m? 0^/2 and Vi{(pi) = m\ The ratio of the masses = nii/m 

is chosen to be smaller than 1 . This choice leads to a two stage inflationary expansion. The first stage is driven mainly 
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by (/), and due to its coupling with the vector field, this stage is anisotropic and it proceeds in a similar fashion as in 
the previous subsection. At the end of the first stage, the field (f>i takes over the expansion, while the field (p starts 
oscillating. Since the expansion of the universe still proceeds, the Hubble friction will drive the field to zero and 
consequently, the coupling f{<j>) approaches unity. Therefore, from we clearly see that the effect of the vector field 
quickly diminishes. Thus, after the first stage, the universe quickly isotropizes and inflation proceeds isotropically 
until the field (f>i starts oscillating. Problems related to an inflationary stage that has increasing anisotropy will be 
eliminated by the two stage modification. 

As we have done in the previous section, we first study the solutions to the background evolution equations, which 
are given in this case as 
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(24) 



where /(0) = e"^"^ /^m^ before. We can obtain slow-roll solutions similar to (fTOl) . in the first stage which proceeds 
anisotropically. Equations (|16p are now modified to become 



6M, 



2 ^/+Ai2 02) 
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3 d 01 



2 2 1 
-TO ^ 01 



(25) 



Initially, the field drives the expansion, so 0^ ^ '/'ij and the slow roll solutions for the first stage is simply given 
similar to (fT9l) as 



2 M2 (c - 

2 TO Adp fj? 



1) 



2 toM„ 



(26) 



and the anisotropy is given precisely as in equation (|2ip. After this first stage ends, the isotropic expansion of the 
universe is driven by 0i and this second stage is described by the following slow-roll solution 



TO Mp /i 



(27) 



We also integrate the system ([M)) numerically and the results are shown in Figs [5]- H) As we have done in the previous 
section, we use the slow-roll solutions and the constraint equation (the first of (j24p ') to determine the initial conditions 
0m I 4>iim ^im ctiri ^ ■ The initial conditions for the scale factors are set as ain — and cri„ chosen to satisfy crit^nd) = 
at the end of inflation. We have chosen 0i„ — 5Mp, 0ii„ = 14Mp and /i = 1/10, which gives around 60 e-folds of 
inflation, with the first 10 e-folds being anisotropic. In this type of background, when the anisotropic expansion 
takes place at the onset of inflation, we expect that only the largest scale perturbations are improved and at small 
scales standard results are recovered. 



III. PERTURBATIONS 



This section is devoted to the study of the perturbations around the anisotropic inflationary background solutions we 
have discussed in the previous section. We will first discuss the decomposition of perturbations around the background 



^ In order to find the approximate slow-roll solutions, one has to ignore the effect of the second field (f>i initially. Although this approxi- 
mation is fairly accurate, we let the numerical evolution to run for a few e-folds until the real attractor solution is reached. 
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FIG. 3: The left panel shows the evolution of the average Hubble expansion. The right panel shows the evolution of anisotropy, 
where H = a and h = a. 




FIG. 4: The phase portrait for the two scalar fields. The straight red line represents (p and the dashed black line represents (jti. 
Time is measured in units of m. 



configuration, which is dilferent from standard scalar- vector-tensor decomposition used for isotropic backgrounds. The 
isotropic background solutions have a 3 dimensional rotational and translational symmetry for constant time slices, 
therefore small fluctuations are decomposed using the representations of the 5*0(3) group. In such a decomposition, 
one identifies decoupled classes of fluctuations and according to their transformation properties, one can identify 
scalar, vector (transverse) and tensor modes(transverse-traceless). Such a decomposition can still be considered in an 
anisotropic background, however fluctuations classified as scalar, vector and tensor would be coupled in the linearized 
calculation ^ (see for example Ell). Instead, we make use of the left over 2 dimensional rotational symmetry 
of the y — z plane of the Bianchi-I metric ^ in order to decompose small fluctuations around the background as 
in [1, [l3|- In such a two dimensional decomposition, there are two decoupled sets of perturbations, the 2d scalar and 
2d vector modes. We will compute the quadratic action of the small fluctuations for both type of modes and study 
them separately. Since the universe isotropizes either after inflation in the single fleld model, or after the first stage 
of anisotropy in the double fleld model, the observables today are best deflned if we use standard isotropic deflnitions 
of fluctuations. Therefore, we will also find the map between the 2d decomposition and the standard 3d isotropic 
decomposition, using which we will determine the power spectra for the gravitational wave polarizations {hx, 
and the curvature perturbation TZ today (which are defined with respect to an isotropic background). This procedure 
is described in detail in the appendix. 



We refer to both the computation of the quadratic action and linearized Einstein equations as the "linearized calculation", since the 
latter follows from the previous. 
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A. Decomposition of Perturbations 



We decompose the perturbations of the metric and the vector field by exploiting the 2d isotropy of the background 



-2$ adix 



abdi {d,B + B, 
P (-2E % - 2 ^^ dj E-d^ Ej - dj E,) j 



(28) 



Additionally, the system of perturbations contain the fluctuations of the scalar field (s) Indices i, j = 2,3 span the 
isotropic 2(i-plane. The perturbations {S(\)cn Xi 5], i?, ag, ai, a} are 2d scalar perturbations, comprising 

one degree of freedom (d.o.f) each, and {-Bi, -Bi, -Ei, OLi\ are 2d vector modes which satisfy the transversality condition 
(9i = • • • = 0). Therefore, the 2d vector perturbations also comprise one d.o.f. each. The 2d vector and scalar 
modes are decoupled at the linearized level. 

We will work in the momentum space where a generic Fourier transformed mode is defined as 



{2nr 



d^kSik) e 



-i kL x—i kT2 y—i fcrs z 



(29) 



We can use the residual 2d-isotropy to fix the comoving momentum to be aligned in the x — y plane without any loss 
of generality. Namely, we set fcrs = and fcr = A:t2 in what follows. Therefore, in Fourier space, the generic 2c? 
decomposition can be written as (before fixing the gauge) 



/ -2"^ -iakLX -ibkrB b B3 \ 

-2a^\E' -abk^kTB -iabk^B^ 

b^ (-2S + 2/c2 £;) ib^kxEs 

\ -252 E J 



6A^{k) = {ao, ai, -ikr a, as) 

Notice that, for the 2d vector modes, fixing kx = kT2 resulted in B2 = B2 = E2 = 0L2 
coordinate transformation + e'^, the metric perturbations transform as 



3 - o(°) e" 



^/iCT ,1/ 



■,(0) ."^ 



where we have decomposed 



e^{k) = {i\-ikLe.-ikTi, e) 



(30) 

0. Under the infinitesimal 
(31) 

(32) 



We choose to fix the coordinate gauge freedom by setting B 
transformation to satisfy 



Y, = E = E^ = Q, namely we solve for the infinitesimal 



a \ a 



r = -E3 



(33) 



Since the action is invariant under the U{1) gauge transformation, we can exploit this symmetry to eliminate one of 
the vector field perturbations. Namely, using the invariance SAfj_ + A, we set a = in ((28| . This choice 

completely fixes the C/(l) gauge freedom. In summary, the complete system of perturbations, after fixing the coordinate 
and C/(l) gauge freedoms, contain three 2d vector modes {B3, a^, B^} and six 2d scalar modes {$, Xi B, ao, ai}. 
We will also have additional scalar modes coming from the fluctuations of the scalar field(s) (50a. We will show in the 
following section that, not all the perturbations are dynamical, namely some perturbations enter into the quadratic 
action without any time derivatives. Such modes are lagrange multipliers and they can be integrated out from the 
action The modes {-B3, 013, ^I*, ai} and perturbations of the scalar fields are dynamical modes. On the other 

hand, {B3, B, ao} are nondynamical modes as we will demonstrate. 

As we have mentioned earlier, it is useful to find the connection between the perturbations defined with respect to 
the 2d decomposition and perturbations defined with respect to standard isotropic decompositions. This is useful in 
discussing the phenomenological consequences. After the universe isotropizes, perturbations evolve in the standard 
way, what is different is that their solutions are nonstandard. However, it is useful to give results in terms of the modes 
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which are commonly used in FRW studies. For definiteness, we find the transformation between our choice of the 
gauge given in equation and the standard longitudinal gauge. Here we present the results and leave the details to 
the appendix. The perturbations in the 2d decomposition are related to the perturbations of the longitudinal gauge 
as 



^3 



k kr 



= 1- 



Ha 



b/ak 



kr 



1 



Ha 

2Hb 



1 

fc2 



b 2kl 



(34) 



In the above equations, the terms in the left hand side are perturbations defined with respect to the 2d decomposition 
and the terms in the right hand side are perturbations in the longitudinal gauge, where , /i+ are the two gravitational 
wave polarizations, "0 is the spatial scalar perturbation and 6(j)^ is the scalar field perturbation (any generic scalar field 
in the model). These relations are useful in computing the power spectra of the gravitational waves and curvature 
perturbation. 



B. General Properties of Coupled System of Perturbations 

In this subsection, we discuss the general properties of the coupled system of perturbations that arise in this study. 
Our main discussion is based on the quadratic action for perturbations, which we need in order to determine initial 
conditions by canonical quantization. We compute the action for perturbations in quadratic order by inserting the 
decompositions (f28|) . 4> = (j>o+5(j) and ipa = 0ao +^0a into ll]), after fixing the U{1) and coordinate gauges as described 
in the previous subsection. We note that the 2c? scalar and 2d vector perturbations decouple in the quadratic action. 
We will provide explicit expressions for the scalar and vector actions in the following sections. Next we integrate out 
the nondynamical modes from the action, to obtain an action in terms of the dynamical modes only. More specifically, 
we extremize the quadratic action with respect to the nondynamical modes and obtain equations of motions for them 
(which are actually constraint equations). Then we solve for the nondynamical modes from their equations of motion, 
and insert the solutions back into the starting action. The final action we obtain depends only on the dynamical 
modes. The next step we perform is to find linear combinations of modes that canonically normalize the quadratic 
action. More specifically, the quadratic action after integrating out the nondynamical modes is formally of the type 
S J dtd^k (Si Kij 6j + Si Ay Sj — Si flfj Sj), where Si denotes the dynamical modes. We find linear combinations 

Si — Rij Scj such that Si Kij Sj = ^ SdSd + . ■ ■ where the dots denote mixed terms. The modes 5^ are the canonical 
modes. At the moment we will study a generic action obtained after canonically normalizing the starting action, 
which will be relevant for both 2d vector and scalar modes. This action reads 

Scan^^ J dtd^fc 1*^*+ (i*^X* + h.c.) -p2$t 1^1 + ^^ $1 (35) 

where # is a column vector made from the canonically normalized perturbations, X is a symmetric matrix and uj'^ 
is a hermitian matrix. In the above action, p is the physical momentum with ~ p\ + Pt where p^ — k^ja and 
Pt = kx/b. We would like to perform a unitary transformation in the field space in order to get rid of the mixed 
terms X ^ and reduce the action to a canonical form that can be used to quantize the fields (and in turn obtain 
the initial conditions) . This can be achieved by a unitary matrix U so that = U '4' (The mixed terms can also be 
eliminated in the level of the Hamiltonian by canonical transformations. See for example js^). We choose the matrix 
U to satisfy 

U + iXU = 0. (36) 
Inserting the transformation back into ([55)) and using (with the unitarity condition on U), we get 



S=\j dtd^k *t*-*^Ut (p2+t^2^xtX) U* . (37) 



We will show in the following sections that the matrix X is of the order of the average Hubble scale H (i.e. X = 
0(H)) and the matrix lo^ is of the order of Hubble scale squared (i.e. = 0{H^)). This means that in the deep 
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UV/subhorizon regime, where p ^ H, these matrices will become negligible compared to p^. Therefore, for modes 
that are deep inside the horizon initially, the above action reduces to a form which has diagonal mass terms 



S =1 I dt(fk 



2 , -p2^t^ + o(7J2) (38) 

In this limit when the mass terms are diagonal, the mode expansion for field ^ simply reads 

The creation and annihilation operators a(fc), a){k) obey the standard bosonic commutation relations and the mode 
functions v(t) satisfy 

i) + u = 0. (40) 
The above equation can be solved using the WKB approximation in the adiabatic limit, which gives 



i;(t) = -i=e-'/-'^*' (41) 
v2a; 



where lo is given by 



to 

r2 



4p2 8 p3 



= p + 0{H^) (42) 

In the above equation, . . . represent higher order corrections which are suppressed by higher powers of H in the 
adiabatic solution. This is also consistent with the accuracy of the mode expansion ([39| . It is of course possible to 
neglect the mixed terms proportional to X and from the beginning. In this way the action for the original fields 
$ is canonical, however we would be working with an accuracy 0{Hin/pin) rather than 0{Hf^/pf^). This situation 
is more difficult in terms of the numerical computation needed, since to achieve the desired accuracy, one has to 
start the numerical evolution much deeper inside the horizon. Moreover, initial conditions in the standard isotropic 
case are also given to an accuracy of 0{Hin/Pin)^ , therefore we choose to perform the field rotation to eliminate the 
X ^ 0{H) terms and work with the fields ^. The field rotation can be performed trivially for the case of 2d vector 
perturbations. The case for the 2d scalar modes are more subtle, therefore we leave the discussion of scalar modes to 
a separate publication. 

The canonical quantization procedure outlined above determines the adiabatic initial conditions for the modes 
^. We are however interested in the values of the entries of $ evaluated at the end of inflation, which are related 
to the primordial power spectra. We will show in the next section that the matrix X along with the off-diagonal 
entries of is proportional to the anisotropy, so they vanish before the end of inflation in the two field model. 
Therefore, it is possible to choose the transformation matrix U that satisfies U(tend) = 1, which in turn implies that 
(w^ + X^ X) U is also diagonal at tend {tend is the time when inflation ends). Then we have $(tend) = '^{tend) 
and we use the following equality to compute the two point functions for the starting modes 

Ite„.=<'i'l*fc> (43) 

for the components of the fields. Therefore, instead of working with the original fields $ we work with the transformed 
fields ^. Notice that passing to the fields ^ was an important step, since the action is in its canonical form with 
diagonal mass terms p^ (up to accuracy 0{H^)) in the deep UV regime. Therefore, this field can be quantized and 
p represents the energy of the each particle created by a^{k). The quantization procedure also determines the mode 
functions v, which in turn determines initial conditions for the perturbations. 

In the next subsection, we will discuss the evolution of the 2d vector modes, which will be related to the power 
spectrum of the gravitational wave polarization hx- 



C. 2d Vector Perturbations 



We now discuss the 2d vector perturbations of the two field inflationary background solution discussed in the 
previous section. As we have argued earlier, the quadratic action for the 2c? vector and scalar modes decouple, so we 
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concentrate on the vector piece here. The action for the 2d vector perturbations in momentum space, up to a total 
time derivative, is calculated as 



S2dV 
^2dV 



2e 



M2 



"3 



b; B3 + h.c. 



2 e -^TTTT- U Q^a -D3 + h.c. 



Af2 

p 



Pi PT-9a' 



Pi 



B, 



3e 



-^pI ( 



Bl B3 + h.c. 



(44) 



where we have used the physical momenta pl = k^/a, pT = kr /b and p = yjp\ + Px- can be seen from (1441) . the 
mode B3 is nondynamical, so it can be integrated out from the action by solving its equation of motion. The solution 
for i?3 is obtained as 



B3 = e' 



Q-2<T 



Pl 



2 ^ PL PA „-2a+CT „, I ^i. / R OAR 



We insert the solution (I45p back into the action and up to a total time derivative, we obtain 



(45) 



72dV 



dtd^kl IffvP + |A_P 



PA PT 



( 



V2AfpP /(./)) 



imA-+iH^ A*_+h.i 



where the canonical modes Hx and A_ are related to the starting modes as 

V2p 



S3 



PL PT Mp 



■ e 2 



a+2a 



g-a/2+(T 



and the mass terms Q,^ are defined as 



nl, = p2 _ 9 2 3 

" ^ 4 4Af2 



!.2 + ^2 



rj2 
"12 



■^12 



ipAPT ( -iPr -^Pl ■ , • 



9 /^Pj , . pIpt Pt 
2 

/'(0) 



Pi 



2M2/2(0) p' 



Pi Pt 



o2 2 1 • 2 I 1 

= P ~4" +4M| 



V2Mppf{(b) V P- 
2 I 12 



p'i 5pI+p't Pa fi't^f 



2{a + a) </' + F'(0) 



/(0) ^ 



2Af 2/2(0) p2 /2(0) J2(0) 



(46) 



(47) 



(48) 



The action (|46p has regular kinetic terms, and therefore there is no ghost instability in the 2d vector sector, as 
expected. Note that the matrix X along with the mass matrix terms 5^2 satisfy 



^1, = 0{H^) , f72,^fi2,=p2 + o(i/2) , 

PA 



= J6d2 - 6a2 _ _L (^(0) + ^^(0^)) „ _L (^2 + 02 I ^ o(i/) 



Mp/(^) 

therefore, the action (j46]) is formally of the type given in ([35|) with and matrix X given by 

PA Pt 



A_ 



, X = 



A 
A 



A = 



V2Mp/(0)p 



(49) 



(50) 
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We can therefore determine the transformation matrix U, that ehminates the mixed terms X $ by solving U + 
I X U = 0, which can be solved by diagonalizing the matrix X. Notice that X R = Dx where 

R.i,f;-M , Z,,.fi _M (51, 



Therefore ([M]) is solved by 

U = Re"*^' ■"■'^'**'Wo (52) 

where Wg is a constant unitary matrix, which does not have any physical relevance. Therefore, we choose Wq = R^, 
so that U = 1 when A = (i.e. when anisotropy vanishes). We can now determine the new mass matrix for the 
transformed fields ^ , which are given by 

nl = [VL^ + X^ X] U (53) 
where the entries of are explicitly given by 

n = A2 + ^^l±3l + cos [21{t)] + CjI^ sin [21{t)] 

^1^2 = >? + ^Hlii - cos mt)] - u^l, sin [2Iit)] 

nl 12 = -nl 2i=i < cos [21{t)] " ^ sin [21{t)] ^ (54) 

where 

/■',,,, I constant t > tiso 

Here, ti^o is the time when the universe isotropizes so X{t > tiso) = 0. The lower limit in the integral is chosen 
for each comoving momentum k such that the constant value for t > tiso is unity ^. This ensures that XJ{tiso) = 
V{tend) = 1 since 

/ cos[X(t)] -isin[J(t)]\ 
- l^-zsin[I(i)] cos[J(t)] j ^^^^ 

Finally, the fields ^ satisfy the evolution equations 

^-a + rf^fc^fc-O (57) 

with initial conditions on the entries of set by the form of the mode function v{t) which are given deep inside the 
horizon as 



1 



2p 



^ ^ -ip-I-\+0{H^) (58) 



V2p V 2p 
up to a nonphysical phase. 

In the next section, we will numerically evaluate the evolution equations (j57p with adiabatic initial conditions (j58p 
for a range of comoving momenta, and determine the two point functions < ^'j^ ^'t > \t^„d which coincides with 
< <^\(^b > |te„d we argued before. Finally, using (|47)) and (pi]) we deduce the power spectrum (which is related 
to the calculated two point functions) for the gravitational wave polarization h^- We will provide the study of the 
2d scalar perturbations in the appendix and leave the study of the power spectrum for the scalar modes and the 
comparison with observations in a separate publication. 



It is possible to keep the lower limit in the integral arbitrary, since it has no physical relevance. However, for an arbitrary value, the 
mass terms for the fields ^ will in general be nondiagonal at the end of inflation. In this case, one has to further rotate the fields (by 
a constant matrix) to obtain the physical modes. We perform this stop from the beginning by choosing tj. corrrectly, so that U = 1 at 
the end of inflation. 
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IV. POWER SPECTRUM 



This section is devoted to the study of the power spectrum of the gravitational wave polarization hx- We will use 
the results of the previous sections to compute the evolution of the fields using their equations of motion given 
in (|57p with the elements of the mass matrices given in (|54p . We perform the integration for a range of co-moving 
momenta and fcy, which we parameterize as 



kL = ^k , kr^y/l-ek , k=^kl + k^ , 0<C<1. (59) 

Therefore, ^ is the cosine of the angle between the comoving momentum and the longitudinal component fc^. For 
each comoving mometa labelled by the two numbers {fc,^}, initial conditions are given in equation (|58p . The initial 
time tin is chosen to be the moment when 

^ 50H{tin) or fc = 50a(i,„)e"(**") (60) 



a{tin) 

so that the physical momentum p is sufficiently larger than the Hubble scale and the adiabatic initial conditions 
(j58p are accurate. Since the mass matrix depends on the background quantities, we integrate the equations 
(|57p together with the background equations (|24p with initial conditions and parameters mi, m2 set as explained in 
subsection B|) . Using p4|l and (j47p we can relate the gravitational wave polarization hx to the canonical field Hx 
at the end of inflation (when cr = a = 0, so a = b) 

% Mp 

We define the power spectrum for the gravitational wave mode hx SiS 

Pxik)^{k^\hx\^)\t=t.„, (62) 

At the end of inflation, the fields ^' is identical to the starting fields so the power spectrum of the gravitational 
wave mode hxis obtained by evaluating the correlation function < '^l >=< $1 > and using (pT|) which gives 

— > 5>P 

Px(fc) = ^e~3"(*=-)|vl/i(ie„.)p. (63) 

In Figures I5I6I we show the power spectrum obtained from the numerical integration of the background and perturba- 
tion equations. In Fig [5] , parts of the power spectrum is shown for fixed values of ^ = 0, 0.5, 1 and 10~^ kiso ^ k < 
10'^/^ kiso- We have deflne kiso to be the value of the comoving momentum, when k ~ a{tiso) H[tiso) = a{tiso) e"^*""-* 
where tiso is the time of isotropization (we define it to be the time when h/H = 10^^ numerically). In Fig. [SI we 

show contour plot of the full spectrum Px{ k). At large scales for k > kiso-, there is more power and the spectrum 
has angular dependence. This in turn could be related to the CMB anomalies. Note that the spectrum reduces to the 
standard nearly scale invariant form when k > kiso- Indeed, we have checked that for k > kiso, the power spectrum 
reduces to 

^x(fc>fc^.o)=^Tfc"- , nT = -2e = -^^ (64) 
where e is the standard slow-roll parameter for the field 0i. Indeed, the the spectrum can be written as 

^ ^ ^iso 

The form of the function /(^, k) is unknown, but it is possible to consider a generic expansion of the form 




LogioPx = LogiQflT + (ci +diC')x + (C2 + d2 C) + (ca + ^3 + 

k. 



Logio ( A ) (66) 
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f=0. f=0.5. f=l 




-1.0 -0.5 0.0 0,5 1.0 1.5 



LogioCk/Mso) 

FIG. 5: Parts of the power spectrum for hx are shown for ^ = (straight line), ^ — 0.5 (dashed line) and ^ — l(dotted line). 

/^x(k) 




-1.0 -0.5 0.0 0.5 1.0 1.5 



Logio[k//:iso] 

FIG. 6: The contour plot of the full power spectrum Px (k). 

We fix the unknown parameters Logj^QOT (which coincides with the Logj^pOT of the spectrum for k < kiso) and 
Ci, C2, C3 by fitting to the numerical spectrum when ^ = 0. Then, we fix the parameters di, d2, c?3 by fitting to the 
numerical spectrum when ^ = 1 (using the previous values of Logj^gOT and ci, C2, C3). Finally, we obtain the values 
of Til, n2, ns by fitting to a spectrum at an intermediate value of ^. Increasing the number of terms in the expansion 
(|66p . increases the accuracy of the fit. For the spectra we have obtained above, the unknown parameters are obtained 
as 

LogioaT ^ -0.43, a ~ -0.23, C2 ~ 0.83, cg ~ 0.38 
di ~ 0.82, da ~ 1.21, dg ~ 0.69 

ni ~ 3.94, 712 ~ 3.38, 713 ~ 3.34 (67) 

In Fig. [71 we show the fitted and numerical power spectra for k < fciso, which are in good agreement. In Fig. [51 we 
show the numerical spectra for k > kiso , which is of the standard nearly scale invariant form. The power spectrum 
for the gravitational wave mode /ix is not of the form Px — Px,iso(k) + (5Px(fc,^) as one would naively expect for 
small anisotropy for k < kiso, with \SPx \ «C |Px,iso|- This is due to the fact the the anisotropic background attractor 
is not continuously connected to standard isotropic solutions, i.e. there is no parameter with a certain limit which 
takes the anisotropic solution to the isotropic one. This is in contrast for example to the ACW model [l5|, where in 
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FIG. 7: The left panel shows the fitted spectrum using (|66|l and the right panel the numerical spectrum for k < kiso- The 
fitted and numerical spectra are in good agreement. 




FIG. 8: The power spectra at small scales for k > kiso, which is of the standard nearly scale invariant form Px <^ fc"^. 

the limit {nri? is the norm of the VEV of the vector field in the model) , the background solution reduces to 

an isotropic hmit. For the current case, attractors with c > 1 (anisotropic) and c < 1 (isotropic) are disconnected in 
the hmit c ^ 1. Therefore, the behavior of the perturbations are dramatically different for k < kiso and the spectrum 
is largely enhanced in this region {\dPx \ ^ |f'x,iso| is not satisfied). For phenomenological considerations, this result 
is not a desired one, since the CMB data would be consistent with a statistical anisotropy at large scales that would 
satisfy \SPx \ ^ l^^x. isol- However, this model is the only vector field driven anisotropic inflationary model that does 
not suffer from any instability to our knowledge. Therefore, we believe that the results obtained in this section are 
still interesting. 

V. CONCLUSIONS 

In this work we have considered anisotropic inflationary models with coupled vector and scalar fields. The coupling 
of the scalar field to the vector field is chosen to preserve the U{1) gauge symmetry of the Lagrangian, therefore, 
the longitudinal polarization of the vector field, which has been shown to cause instabilities [ij, ISOli does not 
exist. The anisotropic expansion is achieved when the vector field has a VEV along one of the spatial directions. The 
anisotropic solution is an attractor, with anisotropy increasing towards the end of inflation, so these types of models 
are counter examples to the cosmic no hair conjecture |35i] . We have considered a generalization of the original model 
introduced in reference [l^], to include additional scalar fields, and specifically concentrated on the double scalar 
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field case. This modification was done in order to achieve isotropization of the universe before the end of inflation. 
The original model with a single scalar field coupled to the vector field has an anisotropic inflationary background 
which persists until the end of inflation. In this type of background, perturbations at all scales are modified. In the 
two field modification, the extra field is not coupled to the vector field, and it is responsible for the overall isotropic 
expansion of the universe only. We have chosen the ratio of the masses of the scalar fields such that the inflationary 
expansion takes place in two steps, the first being anisotropic and the second is isotropic. Therefore, only the largest 
scale perturbations are modified and at small scales the standard spectrum is recovered. As we have discussed in 
the introduction, vector fields with a nonvanishing VEV have been introduced since isotropization takes place very 
quickly in a Bianchi-I background, (when anisotropy is only coming from initial conditions) and thus leads to a fine 
tuning. We have introduced a two phase inflation, with the second phase being isotropic. The transition to the 
second isotropic phase is indeed very quick, however the fine tuning has been relaxed compared to the previous case. 
The second isotropic phase should still be tuned to last around 60 e- folds, but this tuning is not so strict, since 
it only affects the scale where the power spectrum becomes statistically isotropic. Previously, the tuning was more 
stringent, since inflation must last for 60 e-folds in order that the largest scales are modified from the initial anisotropic 
conditions. A more important drawback of the inlfationary setup with only initial anisotropic conditions is that the 
initial singularity is very close to the time of isotropization, which is around one Hubble time. Thus, a calculation of 
the power spectrum shows that it becomes divergent at the largest scales @ ®. 

We have studied the perturbations around the background configuration in detail, taking into account all the degrees 
of freedom of the system. Following the formalism developed in references [1, Q , we have classified decoupled sets 
of perturbations around the anisotropic background and studied them separately. Furthermore, we have computed 
the quadratic action for perturbations and showed that the model is stable and consistent. We also found linear 
combinations of the perturbations which canonically normalize the action. This enabled us to quantize the system 
which in turn determined the initial conditions for the perturbations in the adiabatic vacuum. We then numerically 
integrated the equations of motion for the 2d vector perturbations for a range of comoving momenta, which resulted 
in the power spectrum of the gravitational wave polarization h^- We have only computed the quadratic action and 
determined the canonical modes for the 2d scalar perturbations. We leave the study of the power spectra and the 
consequent phenomenology, which in turn is related to the curvature and /i+ gravitational wave spectrum, to a later 
publication. 

The computed spectrum for has certain interesting features, which might have some relevance to the observed 
CMB anomalies. First of all, the spectrum has angular dependence for k < kiso which breaks the statistical isotropy 
at large scales, therefore leads to an alignment of the lowest multipoles as described in [l^. Moreover, the 

spectrum reduces to the standard nearly scale invariant form for k > kiso, so the higher multipoles are not affected 
from the vector field. We have shown that the power spectrum is formally of the form P^^k) = Px,iso{k) +6Px {k,£,). 
We expect that this will be true also for the spectrum of scalar modes. We have also provided a functional form for 
SPx by fitting to the numerical spectrum in equation (|66p . Phenomenologically acceptable modifications to the large 
scale spectrum should satisfy \SP\ <^ \Piso\, however, we have shown that for the model under discussion, this is not 
the case. Although our results are obtained for the gravitational wave mode hx , we expect a similar behavior also for 
the curvature spectrum, which can be compared with observations. At larger scales, the power spectrum is greatly 
enhanced satisfying \SPx \ > \Px,iso\- Although the anisotropy in the background is small (in the level of a percent, 
as can be seen from eguation ipTjl ). the modification to the spectrum at large scales is not small. This is due to the 
fact that the background attractor solution is disconnected from the standard isotropic solution and therefore, the 
behavior of perturbations are dramatically different. There is no parameter in the attractor solution that connects 
it continuously to the isotropic FRW limit. The anisotropy is always small and proportional to 1/c^ and there is no 
finite limit in c which can continuously deform the attractor solution to the isotropic solution. Therefore, the power 
spectrum cannot be approximated as to be the isotropic piece plus a small correction from anisotropy. This issue 
unfortunately limits the phenomenology of the model. 

Even though the difhculties involved, anisotropic backgrounds still have potentially interesting phenomenological 
consequences. Off-diagonal correlations of the aim coefficients of CMB temperature fluctuations arise in such back- 
grounds and can have relevance to the alignment of lowest multipoles. Moreover, at the largest scales, non standard 
scalar to tensor ratio is obtained, which can be tested by upcoming CMB experiments. The gravitational wave modes 
h+ and hx behave differently at large scales, therefore a future detection of gravitational waves might be tested 
against the consequences of an early stage of anisotropic inflation. In this paper, we have provided a solid example of 
anisotropic inflationary calculation of perturbations and computed the spectrum oi hx polarization of gravitational 



This does not mean that the model is unstable, it reflects that nonlinear effects become important, and one has to go beyond linear 
perturbation theory. 
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waves. Although the obtained spectrum has phenomenologically hmited implications, we believe that our results are 
still relevant, since it is the first complete study of power spectrum in an anisotropic inflationary model driven by a 
vector field to our knowledge. In a further publication, we will study the spectrum of the curvature TZ perturbation, 
which could then be compared with observational data. 
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VI. APPENDICES 



A. Appendix I 



In this appendix we derive the relations (|34p. which are used in the main text to compute the power spectrum of the 
gravitational wave mode ■ These relations are also useful when the curvature and h+ power spectra are calculated. 
For this calculation, one needs to analyze the 2d scalar perturbations, which is done in the next appendix. We will 
study the power spectrum for the scalar modes in a separate publication. 

We start with the decomposition of perturbations in the longitudinal gauge, when the universe is isotropic, which 
is given by 



ds^ = - (1 + 20) dt^ + 2a Vi dt dx' + [(1 - % + hij\ dx' dx^ 



(68) 



where i = 1,2,3. ViS are vector modes satisfying diVi — Q and /i^'s are transverse-traceless modes satisfying 
di hij — hii = 0. Since the spatial slices are isotropic, we can fix the comoving momentum to lie entirely along the 

x-direction k = (fc, 0, 0), without any loss of generality. In this case, the longitudinal gauge can be expressed in 
Fourier space as 



'2(j) V2 
-2a^tP 

(-2V + /1+) 



\ 



V3 


i~2^-h+) ) 



(69) 



We can now rotate the momentum from the x-direction, to a; — y plane in order to make the comparison with the 2d 
decomposition. Namely, we perform the following rotation over the comoving momentum: 



kT 



kh 
kT 



—kr 
kL 



(70) 



where k = \/kj^ + k'^. The matrix appearing in (|70p also determines the Jacobian matrix which fixes the coordinate 

transformation rule from the metric with k = (fc,0, 0) (j69p to a metric with k = (k^, kx,0). The resulting decom- 
position will determine the 2(i-decomposed longitudinal gauge. Using the appropriate scale factors a and b where 
necessary to take into account the anisotropy of the universe, the 2d decomposed longitudinal gauge becomes 



/ -20 



SgUk) 



-ab^h. 



(-2V. 



ab!fh^ 



\ 9 {-2^ - h+) J 

where the superscript "L" indicates that longitudinal gauge is used. Comparing (1711) with (|30p . we identify; 



(71) 
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ikkT 



(72) 



20 



We can now find the infinitesimal transformation between our gauge cfioice oi B = T, 
longitudinal gauge using 



E = E:, 



^Jfj,i^ Jfj,L^,(T ^ iJfia ^,1/ .yen/ ^,fi 



,(0) 



(0) 



JO) 



where 



Finally, the infinitesimal transformation between the two gauges are determined by 



and the 
(73) 

(74) 



^ + ^ ) , 6l = - 



ZjA, r\,rp 



b 1 



a k 



b 2B + 



ikkj 



(75) 



Using these transformations, we can relate the power spectra calculated in the 2d decomposition with the standard 
isotropic definitions of power. For instance, using (j73p and (j75p the dynamical modes of the metric perturbations of 
the 2d decomposition are related to the gravitational wave modes h^, hx and -ip (related to curvature perurbation) as 



kkL 

1 - 



fc|, + b/akl 



Ha 



kT 



\ 2Hb fc2 



a-' 



1 + 



b 2fc? 



Moreover, using the infinitesimal transformation 6(f) = S4>^ ^ 4'£,lj ^hc fluctuation of any scalar field in our gauge is 
related to the one in longitudinal gauge as 

1 



B. Appendix II 

In this appendix, we study the 2d scalar perturbations. As we have done for the case of the 2d vector perturbations, 
we insert the metric and vector field decompositions (|28|) and (j) = 0o(O + '^'/'j 4'i = 4'io{t) + S4>i into the action ([1]) 
(with a = 1), after fixing the U{1) and general coordinate invariance as described previously. The quadratic action in 
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momentum space up to a total time derivative is given by 
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(76) 



We can see from the above action that the modes {$, x, B, ao} are nondynamical, and can be integrated out. More 
precisely, we extremize the action (|76p with respect to the nondynamical modes to obtain constraint equations, which 
we solve for the values of the nondynamical modes. We insert the solutions back into the action ([75]) to obtain an 
action purely in terms of the dynamical modes {(50, (50i, ^, ai} up to a total time derivative. This action reads 
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(77) 



where T and Q are given by 
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(78) 



with /(0i) = Exp / c^0i . Notice that the action ([77]) has regular kinetic terms, therefore the 2d scalar 

sector is free from ghost instabilities as in the case of 2d vector modes. The elements of the hermitian matrix $7^ are 
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where we have defined 

V=2p'^a+ {2pl - p"^) a 
The canonically normafizcd modes {Vi, V2, A_|_} arc related to the starting dynamical modes as 
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When the universe becomes isotropic at the end of the first stage, we have ct = (we also rescaled the initial scale 
factors such that cr = when the universe isotropizes). Using the results of the equations (jTlHTSIIM)) . the canonical 
modes in the isotropic limit is related to the perturbations defined in the longitudinal gauge as 
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where v, vi are the standard Mukhanov-Sasaki variables [33| (defined with respect to the action in the cosmic time). 
All the above quantities are gauge invariant. These relations will be useful when computing the curvature (TZ) and 
gravitational wave polarization spectrum. 
Note that the final action ([77)1 is of the generic form ([55]) , with the corresponding matrix X given by 

/ \ 

X= \ g \ (83) 

\:f g J 

In this case, the solution of the equation U + jXU = is more challenging, since the matrix that diagonalizes X is 
given by 

. / g J'/V2 \ /O \ 

R = -T giV2 _gh/2 , i?A- = VWTg^ o (84) 

^■^ +^ V VT^Tg^/V2 Vr^Tg^i^ ) \Q -^r^Tg^ ) 

which is not a constant matrix, unlike the case for the 2d vector modes. However, the time derivatives of R are 
suppressed by slow-roll parameters, therefore an approximate solution for U can be found. We will discuss the details 
and the computation of the power spectra of the curvature perturbation TZ and hj^ in a separate publication. 
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